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1 Introduction 

Geometrothermodynamics (GTD) is a formalism that has been developed during the past few 
years to describe ordinary thermodynamics by using differential geometry [1]. To this end, 
the states of thermodynamic equilibrium are considered as points of an abstract space called 
the equilibrium space £. Furthermore, we associate to £ a Riemannian metric g in which all 
the geometric properties of £ are encoded. In classical thermodynamics, all the properties of 
a system can be derived from the fundamental equation [2]; analogously, it can be shown that 
in GTD the explicit form of the metric g can be derived from the fundamental equation. It is 
then expected that the thermodynamic properties of the system can be represented in terms 
of the geometric properties of £. In particular, the curvature of £ could be associated with 
the internal mechanical interaction between the constituents of the thermodynamic system, 
i.e. the thermodynamic interaction, so that curvature singularities, in turn, correspond to 
phase transitions. 

In the above approach, starting from a particular fundamental equation, GTD provides 
the geometric structure of the corresponding equilibrium space. However, the formalism can 
also be used to generate fundamental equations. Indeed, if the metric g is assumed to define 
an extremal surface embedded in a phase space T (this will be explained in detail in Sec. 2), 
certain differential equations must be satisfied whose solutions turn out to be mathematically 
well-defined fundamental equations. The consequent question is whether this method can be 
used to generate fundamental equations that could be applied to describe physical systems. 
The main goal of this work is to show two particular cases that can be used to construct 
cosmological models in the framework of general relativity. The idea is to derive all the 
thermodynamic properties from the fundamental equations, and to use them as input to 
construct cosmological models. 

Nowadays, the standard paradigm in the late-time description of the Universe is that it is 
homogeneous and isotropic when averaging over large scales, and that today it is dominated 
by two unknown forms of energy: dark energy, which accelerates the Universe, and dark 
matter that clusters by gravitational instability and is responsible for the formation of the 
structures we see at a very wide range of scales in the cosmos. For a review on the current 
status of cosmology see [3] . 



Because of the lack of a fundamental description of these two ingredients, several alter- 
native proposals have appeared in the literature. In fact, the split of the dark sector into dark 
energy and dark matter is arbitrary, because what we measure in gravitational experiments 
is the energy-momentum tensor of the total dark sector, a property that has been called dark 
degeneracy by M. Kunz in [4]; see also [5-11]. In part for this reason, over the last decade 
the models of unified descriptions of the dark sector have played an increasingly important 
role to describe our Universe. The Chaplygin gas [12, 13] and its generalization [14] will be 
of special interest for us in this work. 

In this paper we find first that the different epochs of the standard cosmological model 
can be described in the context of GTD and that these correspond to the simplest case 
of a system with no internal thermodynamic interaction. Thereafter, we consider a second 
GTD system with thermodynamic interaction that turns out to describe a unified dark sector 
fluid which has as a special case the generalized Chaplygin gas. As a bonus, the so-called 
polytropic fluids can be obtained from this GTD fluid in a certain limit. 

This paper is organized as follows. In Sec. 2, we briefly review the fundamentals of GTD. 
Then, in Sec. 3, we present the cosmological model that follows from a GTD system without 
thermodynamic interaction. In Sec. 4, we study the cosmology of the dark sector GTD fluid 
at the homogeneous and isotropic level. In Sec. 5, we work out the linear perturbation theory 
in order to constrain the free parameters of the GTD model with thermodynamic interaction. 
Finally, section 6 is devoted to discussions of our results and suggestions for further research. 
Throughout this paper we use geometric units in which G = c = k B = h = 1, unless otherwise 
stated. 

2 Basic aspects of geometrothermodynamics 

In classical equilibrium thermodynamics [2], the standard method to describe a thermody- 
namic system consists in specifying a set of n extensive variables E a (a = l,...,n), their 
corresponding dual intensive variables I a , and the thermodynamic potential <I>. The integer 
n determines the number of thermodynamic degrees of freedom of the system. For instance, 
in the case of the ideal gas (n = 2), if we choose the internal energy U as the thermodynamic 
potential $, then E a = (S,V) and I a = (T, —P) so that the temperature T is the dual of 
the entropy S and the (negative) pressure P is the dual of the volume V. All the proper- 
ties of the ideal gas are contained in the fundamental equation U = U(S, V) = (e /V) 2 ' 3 
that satisfies the first law of thermodynamics dU = TdS — PdV from which the expressions 
for the temperature and the pressure, i.e., the equations of state can be derived. Since an 
equilibrium state of the ideal gas can be represented by the corresponding values of S and 
V, all possible equilibrium states form a space £ whose points can be represented by the 
coordinates S and V. 

Notice that using the above notation for an arbitrary system with n thermodynamic 
degrees of freedom, the fundamental equation can be written as <£ = <&(E a ), the first law of 
thermodynamics as d& = I a dE a with I a = 5 a bl , and the coordinates of the equilibrium space 
£ are E a . An advantage of this notation is that it can be used with any thermodynamic 
potential and representation. For instance, to write the above example of the ideal gas 
in the entropy representation one only needs to rewrite the first law of thermodynamics as 
dS = (1/T)dll + (P/V)dV so that the thermodynamic variables are now $ = S, E a = (U, V), 
andl a = (l/T,P/T). 



An important property of classical thermodynamics is that it is invariant with respect to 
Legendre transformations, i.e., it does not depend on the choice of thermodynamic potential. 
Indeed, for the description of the ideal gas instead of U one can also use as thermodynamic 
potential the Helmholtz free energy F = U — TS, the enthalpy H = U + PV or the Gibbs 
energy G = U — TS + PV, without changing the properties of the system. The Legendre 
transformations that generate the potentials F and H are called partial transformations 
whereas G is generated by a total transformation. 

The main idea of GTD consists in associating a differential geometric structure to the 
equilibrium space of a given thermodynamic system in such a way that it does not de- 
pend on the choice of the thermodynamic potential, i.e., it is Legendre invariant. To this 
end, it is necessary to introduce an auxiliary structure called the phase space in which 
the equilibrium space is embedded. To be more specific, let us define the phase space 
as the (2n + 1)— dimensional differential manifold T, with coordinates Z = {Q,E a ,I a }, 
A = 0, ...,2n, equipped with the fundamental Gibbs one-form = d$ — I a dE a [15], and a 
metric G that must be invariant with respect to Legendre transformations. The last condi- 
tion is necessary in order to incorporate in GTD the fact that classical thermodynamics is 
Legendre invariant. In this notation, a Legendre transformation is given by 

{Z A } -)■ {Z A } = {$, E a , I a }, (2.1) 

with 

<S> = $-5 kl E k I l , E i = -I\ P = E\ E j = E j , P = P . (2.2) 

Here i, k, I € / and j € J, where I U J is any disjoint decomposition of the set of indices 
{1, . . . , n}. The metric [16] (summation over all repeated indices) 

G=(d$- I a dE a f + A E a I a dE a dI a (2.3) 

where A is a real constant, is the most general metric we have found so far that is invariant 
under partial and total Legendre transformations, and the last term linear in the extensive 
and intensive variables. 

The equilibrium submanifold £ C T is defined by the smooth map <p : £ — > T, or in 
coordinates ip : {E a } h-> {$(E a ), E a , I a (E a )}, under the condition that <p*(0) = 0, i.e., 

d$ = I a dE a , i.e., Ia=g^, (2.4) 

where (p* is the pullback of ip. These equations are equivalent to the first law of thermody- 
namics and the conditions for thermodynamic equilibrium, respectively. We can associate 
with £ the induced metric 



» = « = A l B *J» ftnt ' <2 ' 5) 

in a canonical manner. One of the main objectives of GTD is to find relations between the 
geometric properties of the equilibrium space £ and the thermodynamic properties of the 
system determined by the fundamental equation $ = §{E a ) [2] that, in turn, is specified by 
the map (p. In particular, one expects that the curvature of £ can be used as a measure of the 
thermodynamic interaction. For instance, in the case of vanishing interaction, one expects 
the curvature to be zero. Let us recall that our interpretation of thermodynamic interaction 



is based upon the statistical approach to thermodynamics in which all the properties of the 
system can be derived from the explicit form of the corresponding Hamiltonian [17], and 
the interaction between the particles of the system is described by the potential part of 
the Hamiltonian. Consequently, if the potential vanishes, we say that the system has zero 
thermodynamic interaction and the curvature should vanish. The equivalence between the 
curvature of £ and the thermodynamic interaction has been shown to be true in the case 
of ordinary classical systems, like the ideal gas and the van der Waals gas [16], and black 
hole configurations in different theories (see [18] for a review). Moreover, the curvature 
singularities of £ turn out to correspond to phase transitions of the thermodynamic system. 
The above description of GTD shows that in order to find explicitly the metric g of 
the equilibrium manifold £ one only needs to specify the fundamental equation <3? = &(E a ). 
This means that one needs the fundamental equation to study the corresponding geometry. 
However, the formalism of GTD allows us to generate fundamental equations by using a 
variational principle as follows. Suppose that the equilibrium manifold £ determines an 
extremal surface in T, i.e., the variation of the volume element of £ vanishes: 



y/det(g)cTE = 0. (2.6) 

£ 

Since g is induced by the metric G that depends on Z , it can be shown [16] that this 
variation leads to a system of differential equations 

UZ A = — L= (^t(g)g ab Z A ) + T A BC Z B b Z c J c = (2.7) 

where □ is the d'Alembert operator. Moreover, this variation implies that the thermody- 
namic potential <3? must satisfy a set of differential equations whose solutions can be written 
as functions of the extensive variables <!> = §(E a ), i.e., as fundamental equations. Two 
particularly simple solutions with <J> = S and E a = {U, V} found in [16] are given by 

S = Cl lnU + c 2 lnV, (2.8) 

and 

S = S In ( U 1+a + cV 1+ P) , (2.9) 



where c±, c 2 , a and /? are real constants. 

The question arises whether these functions, which are obtained as solutions of a geo- 
metric problem, can be used as fundamental equations to describe a thermodynamic system 
with realistic physical properties. This question will be treated in the following sections. 

3 The fluids of the standard cosmological model 

The simplest solution with two thermodynamic degrees of freedom (n = 2) is given by 
equation (2.8) In the special case c\ = 3/2 and c 2 = 1, we obtain the Sackur- Tetrode 
equation that is interpreted as the fundamental equation for the ideal gas [2] . This solution is 
the simplest one in the sense that it corresponds to a system with no internal thermodynamic 
interaction. In fact, introducing the eq. (2.8) into the general metric (2.5) with <3? = S and 
E a = {U, V}, we obtain the particular metric 

. / 2 dU 2 2 dV 2 \ , on , 



A straightforward calculation shows that the curvature of this metric vanishes identically, 
showing that the metric is flat. This can be seen explicitly by introducing the coordinates 
d£ = h 1 i 2 cidU /U and drj = A l ' 2 C2dV/V in which the metric takes the Euclidean form 
g = — (d£, 2 + d?] 2 ). As mentioned above, in GTD we interpret the curvature as a measure 
of the thermodynamic interaction so that a flat metric corresponds to the simplest case of a 
system without interaction. 

The first law of thermodynamics (2.4) in the entropy representation can be written as 

dS = ^dU + ^dV. (3.2) 

Then, from the equilibrium conditions (2.4) we obtain the relationships T = U/c\ and P/T = 
C2/V which lead to the equation of state 

P = -P, (3.3) 

c\ 

where p = U/V. To consider this thermodynamic system in general relativity we assume the 
simplest case of a homogeneous and isotropic spacetime that is described by the Friedmann- 
Lemaitre-Robertson- Walker (FLRW) line element 



ds 2 = -dt 2 + a(ty 



^L + r 2(^2 + sin 2 #2 ) 



(3.4) 



Then, if we assume a perfect fluid source with equation of state (3.3), it is clear that the 
different epochs of the Universe evolution can be obtained by choosing the constants appro- 
priately. So, the choice C2/C1 = 1/3 corresponds to the radiation dominated era, C2 = 
describes the matter dominated era, and 02/ c\ = — 1 corresponds to a vacuum dominated 
cosmology. Consequently, the different fluids of the standard model can be described by 
applying the simplest GTD fundamental equation (2.8) in the context of general relativity; 
in other words, the fluids of the standard cosmological model correspond thermodynamically 
to the simplest possible fundamental equation of GTD. 

It is worth noticing that for the fundamental equation (2.8) the heat capacity at constant 
volume is given by Cy = c\. This opens the possibility of considering the dark energy as 
a non-interacting thermodynamic system with negative heat capacity. In fact, for the dark 
energy fluid we obtained that C2/C1 = —1; therefore, we can assume that C2 > which 
results in a negative Cy. Although most physical systems exhibit a positive heat capacity, 
there are systems for which the heat capacity is negative. Among others, these include self- 
gravitating objects such as stars and star clusters [19]. Furthermore, it can be shown [20] 
that systems with negative Cy are never extensive. We conclude that the dark energy fluid 
can be considered as a non-interacting system with non-extensive thermodynamic variables. 
To further investigate this possibility it is necessary to consider non-extensive variables in 
the framework of GTD. We expect to study this problem in the near future. 



4 A unified description for dark matter and dark energy 

In this section we study the fundamental equation (2.9). According to eq.(2.5), this solution 
generates the thermodynamic metric 
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(1 + a) 2 U 2a [acV l+ P - U l+a ]dU 2 



(U 1+a + cV l +Pf 
+ (1 + j3) 2 c 2 V 2 P[f3U l+a - cV 1+(3 ]dV 2 

-(! + «)(! + ^)cU 1+a V 1+l3 [{l + a)U 1+a + (1 + P)cV l+p ]dUdV 



(4.1) 



for the equilibrium manifold 8. The corresponding curvature is, in general, non-vanishing 
and in the particular case a = 1 and /3 = 1 it can be expressed as 

„ 6U A V A c 2 + 4U 6 V 2 c + AU 2 V 6 c 3 + V 8 c A + U 8 tl . 

R = 5 , (4.2) 

s 2 (c 2 v* + mf v 7 

indicating the presence of thermodynamic interaction. In this sense, this thermodynamic 
system represents a generalization of the system with no interaction investigated in the last 
section. The first law of thermodynamics is again (2.4) and the conditions of thermodynamic 
equilibrium lead to 

1_ S (l + a)U a P S c{l + /3)VP 

T U 1 + a + cV l +^ ' T U l+a + cV 1+ P ' { ' ' 

Then, an equation of state can be written as 

P(U,V)=* 1+ ® VP . (4.4) 

We now consider the large scale evolution of a universe filled with standard model 
particles and the dark sector described by GTD; a subindex d shall denote GTD dark, single 
fluid, variables. We write the equation of state of the dark sector (4.4) in terms of the scale 
factor a(t) and its energy density pa 

P d = -Ca-^ a ~^ p - d a (4.5) 

where we used V = Vb(a/ao) 3 . Also, we set the value of the scale factor today equal to one, 
and defined the constant C = — c(l + /3)Vq Q /(l + a). From this equation one can see that 
the specific case a = /3 in the interval < a < 1 corresponds to a (generalized) Chaplygin gas 
[12, 14]. Moreover, if a = f3 = a fluid, often called dark fluid, which gives exactly the same 
phenomenology as the ACDM model is obtained [11, 21, 22], not only at the cosmological 
level, but also at astrophysical scales. This is because the dark fluid which comprises about 
96% of the energy content of the Universe, partially clusters; for details see [4, 11]. 

Polytropic fluids, extensively used in modeling astrophysical objects, are obtained if 
a = /3 in the interval a < 0. We also note that the case a = 1, dubbed variable Chaplygin 
gas, has been studied in the past and has the advantage over the standard Chaplygin that it 
can develop large inhomogeneous perturbations [23, 24]. 

It is interesting to note that the dark fluid model with a = (3 = leads to a thermody- 
namic metric (4.1) whose curvature vanishes identically. This resembles the case of the GTD 
fluid described in section 3 that generates the fluids of the standard cosmological model. 



The continuity equation p' d = —2>H{p d + P d ) (prime denotes differentiation with respect 
to cosmic time, contrary to conformal time, to be used in the next section) can be integrated 
to give 



Pd 



, ,„ iV(i+«) 

i±| Ca -3(«-/3) +CVa -3(l+«) 



(4.6) 



where Cj is an integration constant. It is convenient to recast this expression into the form 

Aa-tic-fi + (1 _ A)a-^ 1+a A m+a) , (4.7) 



P = PdO 

where we defined p d Q as the value of the dark sector energy density today. The constants are 
related by the equations 

C fl + a \V(i+«) 

A = c + c I (i + f3)/(i + a y Pd0 = {i + p c ' 

C= 1 T ^Ap 1 d r, Cj = Pd ^(l-A). (4.8) 

To ensure the reality and positivity of p d at all times, we must impose the condition A > 
that implies the relation c(l + /3)/(l + a) < 0. Notice that for 1 + a < and positive So, 
it follows that dS/dU < and so the possibility of a negative heat capacity arises, as in 
the case analyzed in the previous section. We will not investigate this case in this section. 
Thus, following eq. (2.9), the entropy of the system must diminish as the configuration space 
grows, and as a consequence the GTD dark fluid has a negative pressure which ultimately is 
responsible to accelerate the Universe. 

It is straightforward to calculate the equation of state parameter of the GTD dark sector 
fluid (w d = Pd/pd), giving 

which has the following behavior 

w d (a -► 0) — ► 0, 

w d (a ^oo) — ► -i±A (4.10) 

1 + a 

w d (a = l) = ~\^A. (4.11) 

1 + a 

Figure 1 shows the evolution of w d as a function of the redshift z = 1/a — 1 for different 
combinations of a and (3; A is kept fixed to the value A = 1/(1 + ^dm/^a) — 0.76, with 
Qi = 8ir Pi0 /3Hl 

Now, the Friedmann equation is given by 

H 2 = —(p d + Pb + p 7 ), (4.12) 

where H = a' /a is the Hubble factor. The energy densities of baryons (pb) and relativistic 
components (p y ) redshift as pb = /Ofeo a_3 and p 7 = p y oa~ 4 , respectively. 




Figure 1. Evolution of Wd as a function of the redshift z. The solid (black) curve corresponds to 
a = /3 = (the ACDM model); the dashed (blue) curve to a = 0.06 and /? = -0.06; the dotted 
(red) curve to a = 0.1 and /3 = 0.2; the dash-dotted (gray) curve to a = /? = 0.5 (a Chaplygin gas). 
A = 0.76 is kept hxed for all the cases. 



To complete with the homogeneous and isotropic description we solve numerically the 
Friedmann equation. We choose the same values for a, (3 and A as in figure 1. The value of 
Pdo is fixed by the flat condition, f^ + ^6 + ^r — 1> giving il^ ~ 0.96. In figure 2 we plot 
the scale factor as a function of the cosmic time for the different chosen combinations of the 
parameter values. 
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Figure 2. Evolution of the scale factor a as a function of the cosmic time t. The solid (black) curve 
corresponds to a = /3 = (the ACDM model); the dashed (blue) curve to a = 0.06 and /3 = —0.06; 
the dotted (red) curve to a = 0.1 and /3 = 0.2; the dash-dotted (gray) curve to a = j3 = 0.5. A = 0.76 
and O c ; = 0.96 are kept fixed for all the cases. The vertical line denotes present time. 



An important quantity for the investigation of the fluid perturbations — to be analyzed 
in the next section — but calculated with purely background quantities is the square of the 
adiabatic speed of sound, c 2 s = Pd/Pd, which can be shown to be 



2 aP d + pp d 
c s = -w d — p— , (4.13) 

or, written as a function of the scale factor, 

2 _ 1 + /3 1 (/3 - «)/(! + a) + /3(1 - .A)fl- 3 ( 1+ ^)/.A 

Cs 1 + a 1 + (1 - ^a-^i+^M (« " P)/0- + a) + {I - A)a~^+P)/A ' l ' ' 

The limits of this expression are c 2 (a — > 0) = 0, and c^(a — > oo) = — (l + /3)/(l + a), ifaj^{3, 
and c^(a — >• oo) = a, if a = /3. This result leads to an important difference between the 
generalized Chaplygin model and the extension found here with GTD. At the cosmological 
background level this fact does not have any consequences, but as we shall see, it is of great 
importance when considering perturbations. To not violate causality we require c 2 < 1; 
consequently, further conditions are imposed over the parameters a and (3. 

The particular case of the Chaplygin gas gives c 2 = —aw d , while for the dark fluid, 
c 2 = 0. The assumption that the speed of sound vanishes has been the starting point in 
several works that study the dark fluid model as an alternative to the ACDM [11, 22]. It turns 
out that both models are fundamentally indistinguishable as long as some general conditions 
are imposed beyond the zero order in perturbation theory. Instead of the Chaplygin gas, it is 
possible to consider its natural extension based upon a constant speed of sound, an approach 
adopted in [25, 26]. 

We note that if a < (3, there is a singularity in the speed of sound at a = [(1 + 
a)(l — A)/A{j3 — a)] 1 ' 3 ^"), this coincides with the moment at which the equation of state 
parameter crosses the phantom barrier, w^ = — 1. 

5 The perturbed Universe 

At small scales (nowadays lesser than about 100 Mpc) the homogeneous and isotropic de- 
scription of the Universe outlined in the last section breaks down. In this section we study 
the deviations of the background cosmology up to linear order in perturbation theory. To 
this end, let us consider scalar perturbations in the Conformal Newtonian gauge, with the 
line element given by 

ds 2 = a 2 {r) [ - (1 + 2^)dr 2 + (1 - 2$)5 ij dx i dx j ] , (5.1) 

where r is the conformal time, related to the cosmic time by dt = adr. The matter fields 
perturbation variables are defined through the expressions 



r° = -p(l + 6), (5.2) 

r = -(p + P)v\ (5.3) 

T i j = P((l + n L )8 i j + V i j ), (5.4) 

where II 1 j is the anisotropic stress tensor. The energy density p and the pressure P denote 
background quantities, and are functions of the conformal time only. The vector v l is called 
the peculiar velocity and is related to the four-velocity u^ of the fluid by the relation v % = 
u 1 /vP '. In the Fourier space we define the velocity 9 = —ikiv 1 and the scalar anisotropic stress 
a = 2hnkjW j w /3(1 + w). 



For a general fluid the energy local conservation equations V^T^" = become [27] 



6 = -(1 + w)(9 - 3$) - ZU (-f- - u 



and 



= -H(l - 3w)6 



w 



± S -^l k *6 + kH-k\ 



(5.5) 



(5.6) 



1 + w 1 + w 

where 5P = Pttl, Sp = p6, H = a/ a and a dot means derivative with respect to conformal 
time. Note that the adiabatic speed of sound can be expressed as c 2 = w — w/3(l + w). To go 
further on, we make the assumption of a perfect fluid, obtaining no anisotropic stresses, a = 0, 
so that the gravitational potentials coincide, <£ = 1 J/. Moreover, if we consider only adiabatic 
perturbations, then the (gauge invariant) entropy perturbation is zero, F = tt^ — c^S/w = 0, 
and the equations for the GTD dark sector fluid become 



5 d = -(1 + w d ){9 d - 36) - m (c 2 s - w d ) 5 d 



and 



-H(l-3c 2 s )e d + 



cjk z 5 d 
l + w d 



+ fc 2 $, 



(5.7) 
(5.8) 



where w d and c 2 are given by eqs. (4.9) and (4.14), respectively. 
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Figure 3. Evolution of the baryonic (left panel) and dark sector fluid (right panel) density contrasts 
as a function of the scale factor a. Solid (black) line corresponds to a — /3 = (the dark fluid 
model). Large-dashed (blue) line to a = /3 = 0.0001. Short-dashed (green) line to a = f3 = 0.0006. 
Dash-Dotted (gray) line to a = 0.0001 and /3 = -0.0001. A = 0.76 and n d = 0.96 are kept fixed for 
all the cases. 



Figure 3 shows the behavior of the baryons and dark sector density contrasts, <5& and 5 d 
respectively, for different chosen parameters a and /3. We note that in the cases with a = (3 
(Chaplygin gases), the density contrasts decay more quickly than those with a ^ /3 (not 
Chaplygin gases). This is because, as shown in figure 4, the squared of the speed of sound of 
the perturbations is positive for the former cases and negative for the latter, enhancing the 
growth of structure. See eq. (4.14) and the discussion thereafter. 

To proceed with the analysis we use the publicly available code CAMB [28] to study 
the anisotropics of the cosmic background radiation. In figure 5 we plot the CMB angular 
power spectrum for different choices of the parameters a and j3, keeping fixed the remaining 
parameters. We note that the larger deviations from the ACDM model show up at large 
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Figure 4. Evolution of the adiabatic speed of sound for the cases considered in figure 4. 



scales. This can easily be understood from the equation of state parameter and the adiabatic 
speed of sound: both of them are nearly zero at high redshifts, thus at early times the GTD 
dark sector fluid behaves essentially as cold dark matter, then at lower redshifts — after 
recombination for the cases shown in figure 5 — they start to diverge from the zero values. 
Consequently, the differences arise mainly through the integrated Sachs- Wolfe effect. This 
enhancement of the low CMB power spectrum multipoles has been found in the past for the 
Chaplygin gas [29], and in general for unified dark models [30]. 
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Figure 5. CMB angular power spectrum for different values of a and f3. A = 0.76 is kept fixed. 

To constrain the parameters of the model, we use the code CosmoMC [31] to per- 
form a Markov Chain Monte Carlo (MCMC) analysis over the eight-parameter space M = 
{Qi,h 2 ,9, T,n s , log A s , a, j3, A}. 9 is defined as 100 times the ratio of the sound horizon to 
the angular diameter distance at recombination, r is the reionization optical depth, n s is 
the spectral index of the primordial scalar perturbations and A s is its amplitude at a pivot 
scale of ko = 0.05 Mpc~ . We take flat priors on the intervals —0.01 < a, (3 < 0.02 and 
0.2 < A < 0.99. 

The observations that we choose to constrain the model are the WMAP seven-years 
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results of the observations of the anisotropics of the CMB [32], and the supernovae type la 
Union 2 data set compilation of the Supernovae Cosmology Project [33]. Moreover, we use 
Hubble Space Telescope (HST) measurements to impose a Gaussian prior on the present 
value of the Hubble constant of Hq = 74 ± 3.6km/s/Mpc [34]. 

Figure 6 shows the marginalized confidence interval in the subspace a — j3; in this 
figure, the region of parameters that corresponds to the Chaplygin gas is represented by a 
solid straight line, and the polytropic case by a dashed line. These lines split the space into 
two regions, a > (3 (with no singular solutions) and a < j3. The circle corresponds to the 
dark fluid (or ACDM) model. 

In figure 7 the 1-dimensional posteriors of the explored space parameter M and the 
derived parameter f^ are shown. For comparison, the results for the dark fluid model are 
also plotted. To translate the latter quantities to the ACDM model language, one only needs 
to use the equations A = 1/(1 + Qdm/&a) and Qd = £Idm + &a, for details see [11]. In table 
1 we present the summary of the results at 0.68 confidence level (c.l.). 

We obtain that the free parameters of the GTD unified fluid have to take values of the 
order of 10 -3 or lesser, although in principle they could be as large as causality allows (for 
the cases a = /3, this is a < 1). This constraints are in agreement with those found for the 
generalized Chaplygin gas in the literature; see e.g. [35, 36]. 




(x10~ 



Figure 6. Contour confidence intervals for the a — (3 subspace of parameters at 68% and 95% c.l. 
The solid line corresponds to the generalized Chaplygin gas, the dashed to a polytropic fluid and the 
circle is the ACDM model. The shading shows the mean likelihood of the samples. 



There is a non-linear effect that we have not considered so far and that arises from the 
fact that in general the relation (P) = P((p),a) does not hold. Therefore, when some scale 
grows and becomes non-linear, the naive averaging procedure is no longer valid. This effect 
has been investigated in the past in [37, 38]; for alternative approaches see [39, 40]. In fact, 
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Figure 7. 1-dimcnsional marginalized probability for the complete set of parameters explored with 
MCMC and the derived parameter $7^. The data used are the WMAP seven-year results, Union 2 
data set supcrnovae compilation and a prior of HST on the Hubble constant. The dashed curves are 
obtained imposing a = ft = 0. 



in our case it follows a relation 

(P) = -Ca- 3 ( Q -« (p)- Q (l - a5 + 2 (<5)), 



(5.9) 



between averages quantities. It is clear that considering these effects complicates the cal- 
culations considerably and it is out of the scope of this work to treat them accurately. We 
expect the free parameters of the GTD fluid to be even more constrained by the corrections 
induced by this non-linear effect. 

6 Concluding remarks 

In this work, we applied the formalism of GTD to construct models of fluids that can be 
used as gravitational sources in the context of relativistic cosmology. First, we considered 
the simplest GTD fluid that corresponds to a thermodynamic system whose equilibrium 
manifold is flat, and found that it can be used to generate all the fluids of the standard 
cosmological model. We also discussed the possibility of considering the dark energy fluid 
as a non-interacting thermodynamic system with negative heat capacity and non-extensive 
thermodynamic variables. 

We then investigated a GTD fluid whose thermodynamic curvature is non-zero in gen- 
eral, indicating the presence of internal thermodynamic interaction. It turned out that this 
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Table 1. Summary of constraints. The upper panel contains the parameter spaces explored with 
MCMC. The bottom panel contains derived parameters. The data used are the WMAP seven-year 
data, Union 2 compilation and HST. 



Parameter 


Best fit a 


io 2 n b h 2 


9 9^1 +0-163 

z.zoi_ 123 





i.038±s : gg; 


T 


0.0892±g : °°» 


10 3 a 


S-ZltlM 2 


10 3 /3 


-1.56±l;li 


A 


0.768+0-™ 


n s 


0.963±°;^ 


logllO 10 ^] 


3.075+g;gg 


ttrf 


Q55+o cos 


to" 


13 84+ - 23 


H a c 


70 41+ 5 - 84 

1 U.^-L_ 3 79 



a. The maximum likelihood of the sample. The quoted errors show the 0.68 c.l. 

b. The Age of the Universe (to) is given in gigayears. 

c. Ho is given in Km/s/Mpc. 



fluid leads to a new cosmological model whose equation of state contains as special cases 
the generalized Chaplygin gas, the dark fluid model, and the polytropic fluids. We showed 
that it is possible to interpret this new GTD fluid as corresponding to a unified model for 
dark matter and dark energy. To prove this, we used the Friedmann equations to perform 
a detailed analysis of the behavior of the state parameter of the GTD dark sector and of 
the corresponding scale factor. The obtained results are in accordance with current cosmo- 
logical observations. The main difference between the generalized Chaplygin gas and the 
GTD fluid consists in the behavior of the adiabatic speed of sound. Although at the cos- 
mological background level this difference does not lead to any considerable consequences, 
the perturbation of the background cosmology shows an essential difference at the level of 
the density contrasts. The square of the adiabatic speed of sound is always positive for the 
Chaplygin gas model but negative in general for the GTD fluid, leading to an enhancement 
of the structure growth in the latter case. Moreover, the analysis of the CMB angular power 
spectrum shows that deviations from the ACDM model appear only at large scales. Finally, 
we find the best fit parameters of the GTD fluid by using current observational data and 
show that the parameters a and j3 must be of the order of 10~ 3 or lesser. 

We conclude that from GTD it is possible to obtain fundamental equations for ther- 
modynamic systems that can be used to develop physically reasonable cosmological models. 
In this work, we analyzed only two simple GTD fluids. It would be interesting to study 
more complicated GTD solutions and their interpretation in the framework of relativistic 
cosmology. 

The microscopic nature of the GTD dark fluid is unknown, as much as the dark matter 
and dark energy in the ACDM model. In this work we have focused on its geometrical 
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description by using the formalism of GTD. 
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